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We investigate the time evolution and stationary states of a stochastic, spatially discrete, pop- 
ulation model (contact process) with spatial heterogeneity and imposed drift (wind) in one- and 
two-dimensions. We consider in particular a situation in which space is divided into two regions: an 
oasis and a desert (low and high death rates). Carrying out computer simulations we find that the 
population in the (quasi) stationary state will be zero, localized, or delocalized, depending on the 
values of the drift and other parameters. The phase diagram is similar to that obtained by Nelson 
and coworkers from a deterministic, spatially continuous model of a bacterial population undergoing 
convection in a heterogeneous medium. 
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I. INTRODUCTION 

In a series of works David Nelson and collaborators 
Q, investigated the extinction versus survival of bac- 
terial colonies undergoing convection in a medium with 
heterogeneous distribution of nutrients. Space is di- 
vided into an "oasis", with high nutrients, and a "desert" 
with low nutrients. The above model is similar to 
models used in mathematical epidemiology in which 
a population is divided into subpopulations (patches) 
and disease-causing factors are non-uniformly distributed 
across patches Q. 

Nelson et al solved the Kolmogorov-Fisher equation 
(KFE) with a convective term, using its similarity to 
a nonlinear Schrodinger equation with a non-hermitian 
perturbation 0, 0, 0] . The eigenvalue spectrum of this 
equation has a diverging localization length with a sharp 
mobility edge. This mobility edge is identified with the 
Fisher velocity, the speed of propagation in the KFE. 
This means that while at high convection velocity all bac- 
teria are transported into the desert and become extinct, 
at low velocity the surviving population is localized near 
the oasis. If the desert can support a non-zero steady 
state population extinction is replaced by a delocalized 
steady state. 

It is tempting to consider the continuum determinis- 
tic equations of Nelson et al Q, Q as a coarse grained 
description of a more microscopic stochastic model of 
this system HSHS@, E3IUE0 A poten- 
tial candidate for such a model is a generalized contact 
process (CP) on a lattice which includes both drift and 
spatial heterogeneity. Using Monte Carl simulations, we 
find that this stochastic model does in fact have a phase 
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diagram which is qualitatively similar to that obtained 
from Nelson's deterministic continuum approach. Var- 
ious transitions between phases (extinct, localized and 
delocalized) are identified. 

Stochastic modeling of the discrete microscopic space- 
time events underlying the type of physical phenom- 
ena considered in 'JJ 0, 0] raises many conceptual and 
methodological questions. These are further complicated 
by the practical limitations of computer simulations — 
the size of the system, the time it can be followed and 
the number of samples are all strictly limited and are 
generally small compared to the macroscopic space-time 
scales of the physical phenomena. Experience neverthe- 
less indicates that judicious microscopic modeling, which 
requires choosing system sizes and simulation times large 
enough to well represent "typical" behavior plus "nor- 
mal" fluctuations, can add to our understanding of the 
macroscopic phenomena 0, IS EH • On the other hand 
the simulation time must not be too large so as to avoid 
rare events such as extinction of populations, which will 
always occur eventually in finite systems. These, how- 
ever, occur on a time scale that grows exponentially with 
the size of the system and so will almost never be seen 
at intermediate time scales. We do not attempt here any 
deep analysis of these important issues but refer to the 
literature HEl [13 • 

The more modest goal of this paper is to investi- 
gate the effect of locomotion with drift, in an inho- 
mogeneous spatial environment, on the behavior of the 
stochastic contact process. This is a frequently used 
model for population dynamics, spread of epidemics, etc. 
[1 El E3- We consider in particular an en- 
vironment consisting of a single oasis surrounded by a 
patch of desert several times the size of the oasis. By 
using periodic boundary conditions, this is equivalent to 
considering a periodic array of such patches which, quite 
unrealistically, are microscopically identical. We first de- 
scribe the basic contact process with additional spatial 
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jumps and find that these jumps have the effect of de- 
stroying correlations. This brings the system closer to its 
mean field behavior. We then investigate the oasis-desert 
system in one and two dimensions. Hydrodynamical scal- 
ing which connects this stochastic model to continuum 
deterministic models is discussed briefly. 



II. DESCRIPTION OF MODEL 

The contact process (CP) is frequently used to model 
population dynamics where thepopulation lives on the 
d-dimensional cubic lattice Z d |E 13 IH IH E3 • At each 
site x of the lattice there is an organism which can exist 
in one of two states {0, 1}: we will refer to as dead or 1 
as alive. Let r\ x = 0, 1 be the variable which specifies the 
state at a site a; of a d-dimensional simple cubic lattice 
Z d and let 77 = {%}, y £ Z d , be the configuration of the 
system at some given time. The stochastic dynamics of 
the process are specified by the following rates: 

i] x = 1 — > rj x = at rate a(x), (1) 
t] x = — > i] x = 1 at rate (X/2d) r/ y (2) 

y:\y—x\=l 

where 2d is the number of neighbors of site x. Equiv- 
alently, a live site dies at rate a(x) or attempts to give 
birth to a new bacterium on a randomly chosen neigh- 
boring site at rate X/2d. It succeeds if that site is empty, 
otherwise nothing happens. We are assuming here for 
simplicity that the birth rate is spatially uniform, al- 
though of course one might expect both to depend on 
local conditions. In particular one expects higher birth 
rates in the oasis than in the desert. What is most impor- 
tant, however, is just the ratio of birth rate to death rate, 
hence the effect we are modeling should be essentially 
captured by just considering variations in the latter. 

As there must be at least one live site to propagate, 
the all-zero state is absorbing, i.e., if r)x_ = for all x at 
any time t it will stay like that forever 0, [n], • When 
a{x) — a is constant, the evolution depends only on the 
ratio A/a. For large A/a a live population can survive. 
As A/a decreases there is a transition from survival to 
extinction at a critical value of this ratio which depends 
on the dimension d of the lattice. In units of a = 1, 
A c ~ 3.30 in one-dimension and A c ~ 1.65 on the square 
lattice HEl 

In the mean field approximation, corresponding to the 
neglect of correlations between different sites, one can 
write a deterministic equation for the density of live sites. 
Thus, for the spatially homogeneous case, with density 
v7s(*)) = p(t), we have 



shows in particular that A 



(MF) 



< A c for d = 1, 2 (or 



dp® 
dt 



-otp + A(l — p)p 2 



(3) 



Starting with some non zero density p(0) < 1 the long 
time solution of Eq. © is p(oo) = if a/A > 1, and 
p{oo) = 1 — a I A, if a/A < 1, and so \i MF ^ = 1. This 



indeed for any finite d). This can be understood by 
noting that correlations between live (occupied) sites are 
positive for the contact process in the stationary state, 
i.e.: (i]xVv) — ( 7 fa)( 7 fe)- Thus correlations inhibit growth, 
since live sites give birth only on unoccupied neighboring 
sites. 

We now add to the CP a mechanism for the move- 
ment of the bacteria via a Kawasaki exchange kinetics 
[1 0, H H ES IH El This permits a bacterium at site x 
to move (jump) to a neighboring site y, \x — y \ = 1. The 
exchanges conserve the total number of individuals in the 
system. These two processes can be considered as follows. 
A clock at each live site x rings at a rate 1, i.e., the prob- 
ability that it will ring in a very small time interval dt 
is equal to dt. When it rings the bacterium at that site 
attempts to jump to a neighboring site with probability 
Pjc, dies with probability (1 — Px)a(x) / c, or attempts 
to give birth to a bacterium on a randomly chosen neigh- 
boring site with probability (1 — Pic)\/(2dc), and with 
probability (l—Pjc)[l — (a(x)+X/2d)/c] nothing happens. 
The constant c is chosen as small as possible to make all 
probabilities non-negative, i.e. c = A/2c? + max^ a(x). 
We note further that if the bacteria decides to jump or 
give birth on a neighboring site which is already occu- 
pied, nothing happens. 

If a jump is decided upon, the appropriate neighboring 
site is chosen as follows: in one dimension a bacterium at 
site x attempts to jump to the site x + 1 with probability 
q and to the site x—1 with probability 1 — q, < q < 1. In 
two dimensions the symmetry of the Kawasaki exchanges 
along the y axis is broken while the exchanges along the 
x axis are symmetric. A bacterium at site (x, y) attempts 
to jump to site (x, y + 1) with probability q/2 and to site 
(x, y — 1) with probability (1 — q)/2 while it tries to jump 
to site (x + l,y) and to site [x — l,y) with equal proba- 
bility 1/4. When q > 1/2 there is a drift to the right or 
up (left or down for q < 1/2). In all cases jumps succeed 
if the target site is empty, otherwise nothing happens. 

When q = 1/2 particles jump in each direction with 
equal probability. This corresponds to a symmetric stir- 
ring. For constant a and q = 1/2 it has been shown 
in Ref. that in the limit, Pk — * 1 one can rescale 
space and time appropriately to obtain rigorously a con- 
tinuum reaction-diffusion equation for the density. A 
space time-scaling, leading to a reaction-diffusion equa- 
tion with drift holds also for q > i; see Section VI. This 
means that when the jumps become very frequent they 
destroy the correlations between sites. Subsequently the 
critical value of A assumes its mean-field value in this 
limit. 

In general the configuration probability distribution, 
p,(rj, t), evolves in time according to the master equation, 



dp-(v, t) 
dt 



P K L K + (1 - P K )L 



G 



l*(v,t) (4) 



where Lq describes the Glauber dynamics, changing rj x 
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to 1 — r] x , while Lk describes the Kawasaki exchanges 
between neighboring sites \x — y\ = 1 [Ell]- 

III. ONE DIMENSIONAL SYSTEM 
A. Homogeneous System 

We carried out Monte Carlo simulations of the contact 
process with a homogeneous death rate, a = 1, and with 
symmetric Kawasaki exchanges for various jump proba- 
bilities Pk- The system was a ring containing 50, 100, 
200, 300, 500 sites (with periodic boundary conditions). 
Results were obtained by taking averages over 10 3 -10 4 
realizations, each starting with an initial configurations 
obtained by putting a bacterium at each site with proba- 
bility pq. For each realization po was chosen uniformly in 
the interval (0,1). As is well known the only steady state 
in the finite system is the absorbing state corresponding 
to r] x — for all x 0, ^| . It is only in an infinite system 
where there can be steady states with a nonvanishing 
population density. It is therefore necessary to extrap- 
olate the critical values of the infinite system from the 
results of the simulations of finite systems p|. The ex- 
trapolated critical values of the CP in one dimension are 
shown in Fig.^as a function of Pk- The simulations in- 
dicate an approximately linear decrease of A c with Pk 
and are consistent with the expectations that stirring 
reduces correlations. As already noted in the limit of 
rapid stirring the system behaves in a mean field way 
with \ C {P K ) ^lasPif^l |. 

We also simulated the homogeneous contact process 
with drift in one dimension. The critical values A c , for 
Pk =0.9 are plotted in Fig. |3 for different values of q. 
We find that the critical value A c decreases approximately 
linearly as q increases to 1; here again the drift helps 
survival by reducing correlations between occupied sites. 
It does this by randomly exposing new empty sites on 
which births can take place. Such an effect does not 
occur for a deterministic ballistic motion which satisfies 
Gallilcan invariance. 



B. Spatial Heterogeneity 

We introduce two types of spatially varying death rates 
in one dimension: 

aj(x) = l + esin(— ) (5) 

/ \ I l-ex6 Oasis ,„s 
aiI ^= { l + cxe Desert ^ 

where x £ {1, 2, N} with periodic boundary condition, 
< e < 1. The oasis is located in {x\l < x < N/5} and 
the desert is everywhere else. For e = 0, a = 1 sets the 
unit of time. 

For both types of spatial variation we used simulations 
to obtain the critical values A c (e), defined as the smallest 



value of A for which a population can survive as a function 
of e. In the absence of exchanges, Pk — 0, the critical 
value A c (e) decreases monotonically as e increases. The 
ratios of A c (e)/A c (0) are plotted against e in Fig. We 
find 

A c (e)/A c (0)~l-e = Mm x a(z)- (7) 

Eq. J7J suggests that the system behaves, for large N, 
as consisting of approximately independent local regions 
(patches) with different values of the death rate a(x). 
If that were true exactly we would have A c (e)/A c (0) = 
1 — e. The stationary population density p(x) is plotted 
in Fig. 2] for both heterogeneous and homogeneous death 
rates using a value of A which is smaller than A c for a = 
1 + e but larger than A c for a = 1 — e. 

C. Spatial Heterogeneity with Exchanges 

In this section we investigate the steady states of the 
contact process with both Kawasaki exchange kinetics 
and a spatial heterogeneity of the form au(x) given in 
Eq. © . The nature of the stationary states of the system 
now depend on the strength of the drift q. 

In the simulations we chose Pk — 0.9 so that we could 
see the change of the stationary states as a response to 
the variation of the drift, while keeping the total simu- 
lation time, which should increase [8j proportionally to 
1/(1 — Pk) manageable. The critical values X c (q) were 
obtained for different values of the drift q in systems of 
size, N = 50, 100, 200, 500 and then were extrapolated to 
those in an infinite system. 

The X-q phase diagram is shown in Fig. [S] for e = 0.5. 
The critical curve (solid line), \ c (q) separates the survival 
and extinction regions. For A < X c (q) the extinct phase, 
i] x = for all x, is the only steady state. For A > X c (q) 
there is a steady state with a finite population density 
which will be approached as t oo, starting with a 
non-zero initial density. There are two survival steady 
states: a localized and a delocalized one. In the localized 
survival phase the population density is non-zero only in 
the oasis and its vicinity while in the delocalized survival 
phase the population density is non-zero everywhere. As 
e decreases the region of the localized phase in the phase 
diagram shrinks. In the limit e = the critical curve 
collapses onto the slanted critical curve of the contact 
process with drift and homogeneous death rate a = 1 
[see Fig. |2. 

The dependence of the stationary density profile p{x) 
on q and A are shown in Figs. and {7\ Fig. shows how 
the population in the oasis is decreased as the wind (drift 
q) increases at A = 1.5. For this value of A the oasis 
is able to sustain a population while the desert is not: 
A C 2(e) < A < A c i(e) where A c i(e) [or A C 2(e)] is defined as 
the critical value for the CP with symmetric Kawasaki 
exchange with Pk — 0.9 and homogeneous death rate 
a = 1 + e [or a = 1 — e]. There is a critical drift, q c (X) 
for each fixed A, above which the population becomes 
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extinct in the long time limit. When the drift is weak, 
[l < Qc (A)] the surviving population is localized near the 
oasis. When the drift is strong, [q > q c (X)] the population 
is "blown" into the desert and becomes extinct. We see 
from Fig. ED (where e = 0.5), that X cl {e = 0.5) ~ 2, 
A c2 (e = 0.5) ~ 0.5 and q c (X c i) ^ 0.6(1). 

Fig. shows the density profile for a strong drift, q = 
0.9, at different A's. The condition for survival, when the 
drift is stronger than the critical drift, is that the desert 
should be able to sustain a population, i.e., A > A c i(e). 
For A > A c i(e) the population survives everywhere and 
its steady state is the delocalized phase. The condition 
for this survival is very much similar to that of the CP 
with homogeneous death rate of a = 1 + e (desert) and 
with strong drift(g = 0.9). 

When q < q c (X), a transition from the localized to 
the delocalized survival phase can be achieved by making 
A cross A c i, the dashed line in A — q phase diagram in 
Fig. [5] When A is very close to A c i , it is also possible to 
have a transition from the localized phase (through the 
extinction phase) to the delocalized phase for constant 
A. This is due to the slanted critical curve in the A — 
q plane: One can draw a horizontal line crossing three 
phases in the A — q plane: the population is localized 
at small drift, (undergoes an extinction transition,) and 
becomes delocalized at strong drift. 



IV. SPATIAL HETEROGENEITY WITH 
EXCHANGES IN TWO-DIMENSIONS 

In d = 2 the exchange kinetics impose a drift along 
the y direction. We consider a periodic box of length 
N and choose the region of oasis to be {(x,y)\2N/5 < 
x < 3N/5,2N/5 < y < 3N/5} and the desert ev- 
erywhere else. We considered square lattices of sizes 
N 2 = 30 2 ,50 2 ,70 2 ,100 2 , with periodic boundary condi- 
tions, for e = 0.5 in Eq. (JSJ and Pk = 0.9, and extrap- 
olated the critical curve to that in an infinite system. 

Note that the steady states in two dimensions also de- 
pends on the shape of oasis. An oasis of narrow width 
along the direction of the drift will yield a very different 
steady states from one having a narrow width perpendic- 
ular to the direction of the drift. One can expect that in 
the former case the steady states will be similar to those 
of the CP with heterogeneous death rate an (x) with 
symmetric exchanges in one dimension while for the lat- 
ter case the steady states will be similar to those of the 
CP with an{x) under the drift in one dimension. 

The X-q phase diagram for the CP with drift and 
with spatial heterogeneity in two dimensions is plotted in 
Fig.[S]for various system sizes with the oasis a square with 
sides of length N/5. The localized [delocalized] phase in 
two dimensions is identified by the presence of a large 
cluster of occupied sites around the oasis, which is "not 
large enough" ["large enough"] to reach both boundaries 
of the system along the x or y direction. The delocal- 
ized, or localized survival phase and the extinction phase 



are separated by the critical curve X c (q). The shape of 
that critical curve is similar to that in one dimension. 
This is not too surprising since, unlike equilibrium sys- 
tems, nonequilibrium systems like the CP, can have phase 
transitions also in one dimension with similar phase dia- 
grams [5J- We can therefore expect similar phase tran- 
sitions between three phases by varying (A, q): examples 
of localized and delocalized phases are depicted in Figs. 03 
andlTTJl 



V. TIME-EVOLUTION OF THE POPULATION 
DENSITY p t (x) 

To study the time evolution of the population density 
Pt{x) undergoing drift in a heterogeneous system, we per- 
formed Monte Carlo simulations in one dimension with 
N = 2000, P K = 0.9 and e = 0.5 in Eq. ©. We selected 
a set of (A, q) parameters from different steady states in 
the X-q phase diagram in Fig. and evolved the system 
from various initial conditions. 

In the localized and in the delocalized phases we 
started with a population density po = 0.01 randomly 
placed. For each set of parameters we simulated 10 4 re- 
alizations of the process with the above initial conditions. 
After 100, 200, 500, 2000 Monte Carlo time-steps we mea- 
sured the population density p t (x) . (In one Monte Carlo 
time step each site is visited once for updating its state, 
on average.) 

Figs.llllandll2lshow the time-evolution of a population 
density under a weak drift (q = 0.59) for a small birth 
rate (A = 1.5) and for a large birth rate (A = 2.5). For the 
small birth rate (A = 1.5) the population grows mostly in 
the vicinity of the oasis and its density profile is skewed to 
the right due to the presence of the drift: under a stronger 
drift the population will die out in this case. For the large 
birth rate (A = 2.5) the population density profile at 2000 
Monte Carlo time-steps decreases gradually in the desert 
as one goes further away from the oasis. Note that in the 
strong drift case, q = 0.9 shown in Fig. the population 
density profile is flatter in the desert for A = 2.4. 



VI. COMPARISON OF THE CONTACT 
PROCESS AND THE KOLMOGOROV-FISHER 
EQUATION 

De Masi et al |8( derived a reaction diffusion equation 
for the symmetric exclusion process with added Glauber 
dynamics by rescaling time and space, x — eja (j — 
—1, 0, 1, ...) and t = e 2 r, a and r being the microscopic 
lattice spacing and r being the time in microscopic units 
(jump time) and then taking the e — > limit with (1 — 
Pk) ~ e 2 . 

We expect that under suitable scalings the model con- 
sidered in this paper leads, when Pk — ► 1, to a continuum 
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macroscopic equation of the form 0, ^| 
dP t (x) 



Of 



= DV 2 P t (x)-vVP t {x){l-P t {x)) 
+ (X-a(x))P t (x)-XP 2 (x) (8) 



where Pt(x) is the bacteria density at time t at position 
x, D is a diffusion constant and w is a convection veloc- 
ity where D and v can be obtained from the microscopic 
model. (Note that on the time and space scales consid- 
ered here there is no contribution to D from the birth 
and death proven, i.e. D — » when Pr — > 0.) 

Eq. © resembles the continuum modified KFE used 
by NelsonQ, @, which generalizes the KFE to account 
for convection and an inhomogeneous growth rate, 



dc(x, t) 
dt 



= D N V 2 c{x, t)-v N Vc(x, t)+U(x)c(x, t)-bc 2 (x, t), 

(9) 

where c(x, t) is the bacteria population density, -Djv is 
the (Nelson) diffusion constant of the bacteria, U(x) is 
the spatially varying growth potential corresponding to 
(A — a(x)) in Eq. JSJ, vn is an externally imposed (Nel- 
son) convection, and b is a parameter that limits the pop- 
ulation number density to a maximum saturation value. 
The only difference between Eq. (JSJ) and Eq. © is the 
nonlinear convection term in Eq. (jHJ). Because our in- 
terest lies mainly in the location of the critical curve, 
we can linearize Eq. ©, while keeping a nonlinear reac- 
tion term. Then the coefficients in the two equations can 
be matched with each other. This makes it possible to 
compare the simulation results with those from the dif- 



ferential equation. We find good qualitative agreement 
between two approaches, even though we don't present 
the results of comparison here. 



VII. CONCLUSIONS 

We studied the steady states and dynamical behav- 
ior of the contact process with asymmetric Kawasaki ex- 
change kinetics and with spatially varying death rate. 
This is a stochastic, spatially discrete, version of the de- 
terministic continuum Kolmogorov-Fisher equation with 
convection velocity, investigated by Nelson and collabora- 
tors 0,0. The phase dia grams from both the stochastic 
and the deterministic cases, birth rate A vs q, are similar: 
two (localized and delocalized) survival phases are sepa- 
rated from the extinction phase by a critical curve, X c (q). 
As is true in many other nonequilibrium situations, the 
difference between one and two dimension is small. Sim- 
ulations of the time-evolution of the population under 
drift in one dimensional heterogeneous systems show how 
a population evolves to the localized or the delocalized 
steady states, depending on the birth rate. 
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FIG. 1: The critical values X c {Pk) as a function of the 
jump probability Pk for q = 0.5. Results are extrapo- 
lated to N = oo (open circles) from system sizes of TV = 
50, 100,200,300,500. The average is taken over 10 3 -10 4 dif- 
ferent realizations starting with random initial configurations 
at density. 
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FIG. 2: The critical value of (q) for a homogeneous system 
as a function of drift q in one dimension. Pk = 0.9 and system 
sizes of N = 100, 200, 250, 500 are used. The function that is 
used for the solid line is 1 — 7(5 — 0.5) with 7 = 0.13. 



0.9 - 

o 

^0.8 
0.7 
0.6 



8\ 


Oa, 




/ a„ 








O 



0.1 



0.2 

e 



0.3 



0.4 



FIG. 3: The ratio of the critical values in heterogeneous sys- 
tem to those in the homogeneous system, A c (e)/A c (0). The 
parameters used for the simulation are Pk = and the re- 
sults are extrapolated to N = oo (open symbols) from N = 
50, 100, 200, 500, 1000. The average is taken over 10 3 - 10 4 



different realizations, 
line is 1 — e. 
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FIG. 4: The stationary density p(x) on an one dimensional 
ring of size N = 2000 for Pk = and A = 4. ai(x) and 
au(x) with e = 0.5 are used, respectively. 




FIG. 5: The X-q phase diagram of the contact process with 
drift and spatial heterogeneity, an (x) with e = 0.5, in one 
dimension. Here Pk = 0.9. The critical curve (solid line with 
filled circles) , obtained by extrapolation to the infinite system 
from N = 50, 100, 200, 500(shown in the inset), separates the 
survival (delocalized or localized) from extinct phase. LP 
stands for the localized phase. 
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FIG. 6: Drift-driven transition from the localized phase to 
the extinct phase. The parameters used for the simulation 
are : N = 2000, A = 1.5, au(x) with e = 0.5 and P K = 0.9. 
Drift varies from q — 0.5 where population is localized in 
the vicinity of the oasis to q = 0.598 where the population 
becomes extinct. 



9 

l i . , . , . , . 1 




FIG. 7: Transition from the delocalized survival phase to the 
extinct phase. The parameters used for the simulation are 
system size of N = 2000, q = 0.9, au(x) with e = 0.5 and 
Pk =0.9. Birth rate varies from A = 2.4 where population 
is extended everywhere to A = 1.9 where the population be- 
comes extinct. 
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FIG. 8: The X-q phase diagram of the contact process with 
drift and spatial heterogeneity in two dimension. The pa- 
rameters used for the simulation are Pk = 0.9, an(x) 
with e = 0.5. The critical curve (solid line with filled cir- 
cles) is obtained by extrapolation to the infinite system from 
N 2 = 30 2 ,50 2 , 70 2 , 100 2 (shown in the inset). The oasis is a 
square of length N/5. LP stands for the localized phase. 
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FIG. 9: The localized phase in two dimension for A = 1.05 
q = 0.5, an(x) with e = 0.5, P K = 0.9 and iV 2 = 100 2 . 
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FIG. 11: Time evolution of population density pt{x) for A = 
1.5, q = 0.59, an(x) with e = 0.5, P K = 0.9 and N = 2000, 
i.e., for a weak drift in the localized phase. 10 4 realizations 
starting with 1 percent of randomly distributed population 
are averaged at given Monte Carlo time-steps. 
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FIG. 12: Time evolution of population density pt(x) for A = 
2.5, q = 0.59, a// (as) with e = 0.5, P K = 0.9 and N = 2000, 
i.e., for a weak drift in the delocalized phase. 10 4 realizations 
starting with 1 percent of randomly distributed population 
are averaged at given Monte Carlo time-steps. 



